Modulation of energy spectrum and control of coherent microwave transmission at single-photon level by longitudinal field in a superconducting quantum circuit
Guo Xueyi1, 2, Deng Hui1, Li Hekang1, 2, Song Pengtao1, 2, Wang Zhan1, 2, Su Luhong1, 2, Li Jie1, Jin Yirong1, †, Zheng Dongning1, 2, ‡
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
University of Chinese Academy of Sciences, Beijing 100049, China

 

† Corresponding author. E-mail: jyr-king@iphy.ac.cn dzheng@iphy.ac.cn

Project supported by the Ministry of Science and Technology of China (Grant Nos. 2014CB921401, 2017YFA0304300, 2014CB921202, and 2016YFA0300601), the National Natural Science Foundation of China (Grant No. 11674376), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB07010300).

Abstract

We study the effect of longitudinally applied field modulation on a two-level system using superconducting quantum circuits. The presence of the modulation results in additional transitions and changes the magnitude of the resonance peak in the energy spectrum of the qubit. In particular, when the amplitude λz and the frequency ωl of the modulation field meet certain conditions, the resonance peak of the qubit disappears. Using this effect, we further demonstrate that the longitudinal field modulation of the Xmon qubit coupled to a one-dimensional transmission line could be used to dynamically control the transmission of single-photon level coherent resonance microwave.

1. Introduction

Superconducting circuits are promising candidates in realizing scalable quantum computing.[1] The quantum coherence properties of superconducting qubits have been improved significantly over the past decade with a five orders of magnitude increase in coherence time, owing to progress in the optimization of materials, device design, and fabrication process.[2] The steady effort on the issues of measurement technique and isolation of modes and infrared radiation has also contributed to the performance of superconducting circuits.

Superconducting qubits can be regarded as two-level systems and the study of their quantum behavior is essential. Quantum two-level systems have been extensively studied in-depth. The simplest situation is its evolution driven by a resonant transverse field. The stationary Hamiltonian of the two-level systems is , where ħω0 is the energy level difference between the ground and excited states and σz is the Pauli matrix. When a weak resonant transverse field drive ħλxcos(ω0t)σx (λxω0) is applied, the rotate-wave approximation (RWA) is valid and the two-level system will do the Rabi oscillation.[3] The RWA ignores the counter rotate part that should not be ignored when the drive is strong.[4] Compared to the situation of driving a two-level system with transverse field, the study on the effect of longitudinal field applied to a two-level system is relatively lacking in the literature. Frequently, one applies a longitudinal bias dc-field to modify the energy level of a qubit. Recently, the effect of simultaneously applying longitudinal time varying field modulation and transverse field drive to a two-level system has been studied both theoretically and experimentally.[515] In this case, the dynamics of the two-level system is not simply described by Rabi oscillation. In a particular case, it is found that if the amplitude and frequency of the longitudinal field modulation (LFM) field take a certain ratio, the Rabi oscillation would disappear due to the modification of the energy levels by the LFM.[10] By using this effect, it has been demonstrated that the coupling strength between a superconducting flux qubit and a co-planer waveguide resonator could be changed in a controlled manner.[16]

In this work, we explore the effect of LFM on a two-level system using a Xmon superconducting qubit device. In Section 2, we consider the situation of simultaneously applying LFM and transverse field drive, and experimentally measure the excitation spectra of the Xmon. The results are in agreement with the theoretical analysis. In Section 3, we investigate the microwave photon transmission behavior of a superconducting one-dimensional coplanar waveguide (CPW) transmission line coupled to a Xmon qubit that is driven by a LFM field. The propagating wave is scattered by the qubit and the transmission is LFM dependent. We demonstrate dynamical control of the transmission characteristic of single-photon level coherent microwave through LFM pulse.

2. Sample fabrication and measurement set-up

Two samples were used in this work. The superconducting qubit is of Xmon type in both samples.[17] The energy gap of the qubit can be adjusted by an external flux bias. In the first one, a Xmon qubit is capacitively coupled to a λ / 4 coplanar waveguide (CPW) resonator that is coupled to a CPW transmission line. In this sample, the qubit state is readout by the dispersive method via the λ / 4 resonator. In the second sample, a Xmon qubit is directly coupled to a transmission line and the interaction between an artificial atom (i.e., the Xmon qubit) and photons transmitted in the transmission line can be investigated. The optical micrographs of the two samples are shown in Figs. 1(a) and 2(a), respectively.

Fig. 1. (color online) (a) Schematic illustration of Rabi oscillation. (b) Schematic illustration for the disappearance of Rabi oscillation with LFM.
Fig. 2. (color online) (a) Optical micrograph of the Xmon qubit of the first sample. The green part is a λ/4 CPW resonator for readout. The red cross part is the capacitor of Xmon and the inset shows the DC-SQUID of Xmon which works as an adjustable non-linear inductance. The XY control line is used to apply a transverse field drive and the Z control line is used to apply the bias field and thus to apply LFM. (b) Schematic diagram of measurement setup. (c) Energy spectra of the Xmon qubit measured as a function of bias voltages. (d) Xmon qubit energy spectra variation with LFM amplitude. The qubit is biased at Φ = 0.15758Φ0, corresponding to an energy difference ω0 = 6.1455 GHz and the LFM frequency is fixed at 20 MHz. (e) Resonate peak intensity as a function of LFM amplitude for 6 branches (corresponding to n = 0, −1, …, −5) shown in (d). The solid lines are fitting of the Bessel function (Eq. (5)) to the experimental data.

The samples were fabricated using a process involving electron-beam-lithography (EBL) and double-angle evaporation. In brief, a 100 nm thick Al layer was firstly deposited on a 10 mm × 10 mm sapphire substrate by means of electron-beam evaporation, followed by EBL and wet etching to produce large structures such as microwave coplanar-waveguide resonators/transmission lines, capacitors of Xmon qubit, and electric leads. The EPL resist used was ZEP520 and the wet etching process was carried out using aluminum etchant type A. In the next step, the Josephson junctions of qubits were fabricated using the double-angle evaporation process. In this step, the under cut structure was created using a PMMA-MMA double layer EBL resist following a process similar to that reported in Ref. [18]. During the evaporation, the bottom electrode was about 30 nm thick while the top electrode was about 100 nm thick with intermediate oxidation.

In the measurements, the sample was mounted in an aluminum alloy sample box which was fixed on the mixing chamber stage of a dilution refrigerator. The temperature of the mixing chamber was below 15 mK during measurements. The input microwave lines and qubit fast bias control lines were heavily attenuated. Lines for qubit dc bias control were filtered using filters (RLC ELECTRONICS F-10-200-R) that functioned as combination of low-pass filter and copper powder filter. A bias-tee was used to combine the dc bias and LFM produced by an arbitrary wave generator (AWG). The microwave output signal from the transmission line was amplified (≈ 39 dB) by a cryogenic HEMT amplifier mounted at the 4 K stage and a room temperature amplifier (≈ 38 dB) before being measured either by a home-built heterodyne acquisition system or a vector network analyzer. The measurement set-ups are schematically shown in Figs. 1(b) and 2(b).

3. Excitation spectra of Xmon with LFM

Consider a simple two-level system subject to simultaneous longitudinal field modulation and transverse field drive. The Hamiltonian can be expressed as

Here, H0 is the stationary Hamiltonian of the qubit, Hlongitudinal(t) is the LFM term, and Htransverse(t) is the transverse field drive. Note that this transverse field drive is not necessarily in resonance with the qubit, i.e., ωe is not necessarily equal to ω0.

To study the effect of LFM, we apply uniform transformation

to Eq. (1). Then the Hamiltonian becomes
where Jn[x] is the n-th Bessel function of the first kind. Here we use the relation
From Eq. (3), we can see that only terms that meet ωe = ω0 + l (n = 0, ±1, ±2, . . . ) will be non-zero over a time average. Thus, in this case, the Hamiltonian is further simplified to
Therefore, we obtain the effective Rabi frequency for transverse field drive at frequency ωe = ω0 + l.

From Eq. (5), we can easily see that the system should show two features. At first, because n can be any arbitrary integer, there should be other transitions besides the original transition at 0. Secondly, equation (5) shows that the effective Rabi frequency depends on the modulation amplitude λz and frequency ωl through the Bessel functions. It is known that the Bessel functions have zeros for certain values of arguments. Therefore, an important feature stemmed from Eq. (5) is that the effective Rabi frequency would become zero for certain ratio of λz / ωl. For instance, in the case of n = 0, while . Figure 1 illustrates the two processes with and without the LFM.

To verify the effective Rabi oscillations induced by the LFM at different amplitude and frequency, we performed excitation spectra measurements on a Xmon qubit by applying the continuous weak transverse field drive through XY control line and LFM through Z control line. Similar experiments were conducted in Ref. [15] for a superconducting transmon qubit. In the spectra measurements, the intensity of the corresponding transitions is proportional to the population that can be obtained from the density matrix.

The time evolution of the density matrix of the system is described by the master equation

where is the Lindblad operator that accounts for the energy relaxation and dephasing process of the system. The final equilibrium population of |1⟩, i.e. ρ11, is related to in Eq. (5), and thus

We first measured the basic characteristic parameters of the samples. For the first sample, the maximum frequency of the Xmon qubit is and its anharmonicity (ω12ω01)/2π = −215 MHz. Here, we label the lowest three levels of Xmon as |0⟩, |1⟩, and |2⟩. The resonator frequency would shift according to the qubit state when they are dispersively coupled to each other. Thus, for the maximum qubit frequency bias point, the resonator frequency becomes ωr,|0⟩/2π = 6.5564 GHz and ωr,|1⟩/2π = 6.5542 GHz when Xmon is at |0⟩ state and |1⟩ state, respectively. The coupling strength between the Xmon and the readout resonator is , where Δω01ωr and ηω12ω01. Finally, by performing time domain measurements, we obtained the qubit energy relaxation time T1 = 5.27 μs.

In Fig. 2(c), we show a part of the energy spectra of the Xmon qubit as a function of dc bias. For the Xmon qubit, the energy difference between the ground state and the first excited state is , where EJ and Ec represent the junction energy and charge energy, respectively. EJ can be adjusted by external flux Φ via the relation EJ(t) = EJ,max|cos(πΦ(t)/Φ0)|. Because ω01 is insensitive to Φ at the maximum frequency point, we bias the qubit to a point Φ = 0.15758Φ0 with dc-bias. The correspoing qubit frequency is ω0 = ω01 = 2π · 6.1455 GHz.

The dependence of the Xmon qubit energy level difference on the external flux is not linear. Therefore, in order to obtain a sinusoidal modulation, we need to calibrate the output voltage waveform of the AWG used in our measurements. This can be easily done with the AWG according to the above mentioned relations among ħω01, EJ, and Φ.

We measured the excitation spectra with a fixed LFM frequency of 20 MHz and investigated its variation with LFM amplitude. The experiment results are shown in Fig. 2(d). Clearly, as would be expected from Eqs. (2) and (3), additional resonate transition peaks appear. The frequency difference between the neighbor peaks is ωl = 20 MHz. It also shows that at certain LMF voltages, the resonate peaks disappear, as would be expected from Eq. (5). In Fig. 2(e), we present the resonate peak intensity as a function of LFM amplitude for 6 branches (corresponding to n = 0, −1, …, −5 ) shown in Fig. 2(d). The solid lines are fitting of the Bessel function (Eq. (5)) to the experimental data. The measured result is in good agreement with the theoretical analysis.

We note that in Fig. 2(d), the excitation peak intensity is weaker in the upper branches than that in the lower ones. We attribute this to the variation of energy relaxation time T1 at different bias points. The data shown in Fig. 2(c) appear to support this explanation with relatively weak excitation peaks at the high frequency region.

4. Control of transmission of single-photon level coherent microwave by LFM

In the previous section, we showed LFM induced transparency with Xmon excitation spectra. In this section, we coupled a Xmon in a one-dimensional transmission line to study the related quantum optic effect when applying LFM to Xmon.

The effect of superconducting artificial atom on the transmission coherent microwave photons has been reported in Refs. [19] and [20]. For example, strong atom–field interaction was realized with a superconducting flux qubit embedded in a one-dimensional transmission line, demonstrating a high degree extinction of propagating wave and showing the resonance fluorescence.[19] Later in another work, a transmon qubit was embedded in a one-dimensional transmission line, and the authors used the ATS effect to control the passage or reflection of single-photon level coherent microwave.[20] Other quantum optical effects and applications due to the interaction between artificial atoms and propagating photons in one-dimension transmission lines were studied theoretically and experimentally in Refs. [21]–[37].

In the following, we show that the LFM on a qubit can change the transmission properties of single-photon level coherent microwave. As shown in Ref. [19], for a qubit placed at x = 0 and an incident microwave V0(x,t) = V0eikx−iω0t, where k is the wavenumber and ω0 is the frequency, if the qubit frequency also is ω0, the scattered wave is Vsc(x,t) = Vsceik|x|−iω0t. Then, the reflected wave is Vr(x,t) = Vsc, and the transmitted wave is the coherent interference of the scattered wave and incident wave Vt(x,t) = V0 + Vsc. Because of the boundary condition for the scattered wave, the scattered wave is in opposite phase with the incident wave. Hence, these two waves are in destruction interference in the transmitted direction.

If the qubit is modulated by LFM, the energy level difference becomes ħ(ω0 + λzcos(ωlt)). It would generate scattered waves with multiple frequency, and the amplitude of each different frequency component varies in the form . This may be understood with the following formula:

Thus, in the transmitted direction, when the incident wave frequency equals to ω0 + l (n = 0, ±1, ±2, . . .), it also will be in destruction interference with the scattered wave, and the net wave amplitude is affected by the scattered wave amplitude. Supposing the frequency of incident wave ωi = ω0 + l and using the opposite phase condition, we obtain the transmission coefficient
In experiment, we can measure how S21 change with the LFM amplitude.

The measured data are shown in Fig. 3. The energy spectra of the Xmon qubit as a function of dc-bias voltage are displayed in Fig. 3(d). To study the effect of LFM, the qubit is tuned at 5.779 GHz while a LFM with a fixed frequency 50 MHz is applied to the qubit via the Z control line. Again, the AWG output waveform is calibrated for producing LFM using the method mentioned in the previous section. The transmission spectra of the sample measured at different frequencies as a function of LFM amplitude are shown in Fig. 3(e). Clearly, the results are in agreement with Eq. (9) and show that the transmission can be modulated by LFM. According to Eq. (9), the transmission coefficient becomes unit when is zero. Thus, we may choose appropriate λz and ωl to obtain full transparency of the incident microwave waves at a frequency in resonant with the qubit.

Fig. 3. (color online) (a) Optical micrograph of the second sample. The green part is a one-dimensional open transmission line. The red part is the Xmon and the inset is the dc-SQUID of Xmon. LFM can be applied through Z control line. (b) Schematic diagram of experiment setup for transmission spectrum measurement. The incident microwave is fed in through port 1 and the output signal is picked up on port 2. LFM is generated by an AWG and the dc bias is generated by a dc voltage source. These two bias signals are combined by a bias-tee and connect to the Z control line. (c) Transmission coefficient S21 as a function of incident microwave power at the frequency in resonant with the Xmon qubit. In this case, the frequency is 6.42793 GHz. The red dot is the experiment data, and the blue solid line is the fitting curve of Eq. (??). (d) Qubit energy spectra variation with dc bias obtained from S21 measurements without LFM. (e) Qubit energy spectra at a fixed dc-bias, corresponding to a qubit frequency 5.779 GHz, measured from S21 measurements with LFM. The LFM frequency is 50 MHz and the amplitude is varied.

In the following, we demonstrate transient control of single-photon level coherent microwave transmission by turning on and off the LFM. At first, we measure the transmittance as a function of the incident wave power and the results are shown in Fig. 3(c). Transmission behaviors of such kinds of samples have been studied in Refs. [19] and [20]. If the incident wave is in resonant with the qubit, its transmission coefficient is

where r0 = η/(1 + 2Γϕ/Γ1). η is a dimensionless coupling efficiency of Xmon and transmission line, and Γ1 is the energy relaxation rate of Xmon from |1⟩ to |0⟩. Γ2 = Γ1/2 + Γϕ is the dephasing rate and Γϕ is the pure dephasing rate. Ω is the Rabi frequency of Xmon and Ω ∝ |V0|. The solid line in Fig. 3(c) represents the fitting of the experimental data to Eq. (10). We define an average number of photons per interaction time 2π/Γ1 as ⟨n⟩ ≡ 2πP/(ħω0Γ1). At a low incident microwave power (≈ −130 dBm and ⟨n⟩ ≈ 2 calculated with Γ1 ≈ 2π· 10 MHz), we measured the transmission coefficient with LFM being turned on and off.

According to the data in Fig. 3(e), we determine that for the 50 MHz LFM, the transparency appears at ω0 = 5.779 GHz when the amplitude of LFM waveform output by the AWG is about 0.027 V. The data are shown in Fig. 4 in which a LFM pulse of 1 μs is applied. Clearly, we observe the modulation of the transmission by LFM. We would like to point out that this is an alternative approach as compared to the previously reported method of using the ATS effect[20] or direct detuning the 0–1 transition frequency. It may find applications in building quantum networks. We note that there is still about 20% transmission when the LFM is turned off. We believe that this could be improved by increasing the coupling strength between the Xmon qubit and the transmission line.

Fig. 4. (color online) Dynamic control of the transmittance of single-photon level coherent microwave with LFM pulse. Incident microwave is at frequency 5.779 GHz. We apply a calibrated 1 μs LFM pulse and measure the response of S21. When the LFM is applied, the microwave is totally transmitted. Without LFM, the microwave is 80% reflected.
5. Conclusion

We have investigated the effect of longitudinal field modulation on the properties of superconducting Xmon qubit. In the first experiment, we measured the excitation spectra of a superconducting Xmon qubit while simultaneously applying a longitudinal field modulation through the Z control line and a transverse field drive through the XY control line. We observed the appearance of transitions at frequencies different from the qubit frequency by multiple times of the modulation frequency. Furthermore, transition peaks disappeared when the amplitude λz and frequency ωl of the modulation field met the condition λz ≈ 2.4 ωl, as if the qubit became transparent to the resonance photons. In the second experiment, we used this effect to dynamically control the transmission of coherent microwave at the single-photon level.

Reference
[1] Devoret M H Schoelkopf R J 2013 Science 339 1169
[2] Oliver William D Welander Paul B 2013 MRS Bulletin 38 816
[3] Rabi I I Ramsey N F Schwinger J 1954 Rev. Mod. Phys. 26 167
[4] Ying Wu Yang Xiaoxue 2007 Phys. Rev. Lett. 98 013601
[5] Z Zheng H 2012 Phys. Rev. 86 023831
[6] Yan Y Z Zheng H 2014 Phys. Rev. 90 053850
[7] Z Yan Y Goan H S Zheng H 2016 Phys. Rev. 93 033803
[8] Yan Y Z Zheng H Zhao Y 2016 Phys. Rev. 93 033812
[9] Yan Y Z Luo J Zheng H 2017 Phys. Rev. 96 033802
[10] Liu Y X Yang C X Sun H C Wang X B 2014 New J. Phys. 16 015031
[11] Silveri M Tuorila J Kemppainen M Thuneberg E 2013 Phys. Rev. 87 134505
[12] Nakamura Y Pashkin Y A Tsai J S 2001 Phys. Rev. Lett. 87 246601
[13] Wilson C M Duty T Persson F Sandberg M Johansson G Delsing P 2007 Phys. Rev. Lett. 98 257003
[14] Wilson C M Johansson G Duty T Persson F Sandberg M Delsing P 2010 Phys. Rev. 81 024520
[15] Li J Silveri M Kumar K Pirkkalainen J M Vepsäläinen A Chien W Tuorila J Sillanpää M Hakonen P Thuneberg E Paraoanu G 2013 Nat. Commun. 4 1420
[16] Wu Y Yang L P Zheng Y Deng H Yan Z Zhao Y Huang K Munro W J Nemoto K Zheng D N Sun C P Liu Y X Zhu X Lu L 2016 arXiv:1605.06747 [quant-ph]
[17] Barends R Kelly J Megrant A Sank D Jeffrey E Chen Y Yin Y Chiaro B Mutus J Neill C O’Malley P Roushan P Wenner J White T C Cleland A N Martinis J M 2013 Phys. Rev. Lett. 111 080502
[18] Kelly J S 2015 Fault-tolerant Superconducting Qubits Ph.D. Dissertation Santa Barbara University of California
[19] Astafiev O Zagoskin A M Abdumalikov A A Pashkin Y A Yamamoto T Inomata K Nakamura Y Tsai J S 2010 Science 327 840
[20] Hoi I C Wilson C M Johansson G Palomaki T Peropadre B Delsing P 2011 Phys. Rev. Lett. 107 073601
[21] Shen J T Fan S 2005 Phys. Rev. Lett. 95 213001
[22] Abdumalikov A A Astafiev O Zagoskin A M Pashkin Y A Nakamura Y Tsai J S 2010 Phys. Rev. Lett. 104 193601
[23] Astafiev O V Abdumalikov Jr A A Zagoskin A M Pashkin Y A Nakamura Y Tsai J S 2010 Phys. Rev. Lett. 104 183603
[24] Abdumalikov A A Astafiev O V Pashkin Y A Nakamura Y Tsai J S 2011 Phys. Rev. Lett. 107 043604
[25] Wilson C M Johansson G Pourkabirian A Simoen M Johansson J R Duty T Nori F Delsing P 2011 Nature 479 376
[26] Hoi I C Palomaki T Lindkvist J Johansson G Delsing P Wilson C M 2012 Phys. Rev. Lett. 108 263601
[27] Hoi I C Kockum A F Palomaki T Stace T M Fan B Tornberg L Sathyamoorthy S R Johansson G Delsing P Wilson C M 2013 Phys. Rev. Lett. 111 053601
[28] Hoi I C Wilson C M Johansson G Lindkvist J Peropadre B Palomaki T Delsing P 2013 New J. Phys. 15 025011
[29] Koshino K Terai H Inomata K Yamamoto T Qiu W Wang Z Nakamura Y 2013 Phys. Rev. Lett. 110 263601
[30] Koshino K Inomata K Yamamoto T Nakamura Y 2013 Phys. Rev. Lett. 111 153601
[31] Loo A F V Fedorov A Lalumi.ere K Sanders B C Blais A Wallraff A 2013 Science 342 1494
[32] Inomata K Koshino K Lin Z R Oliver W D Tsai J S Nakamura Y Yamamoto T 2014 Phys. Rev. Lett. 113 063604
[33] Mlynek J A Abdumalikov A A Eichler C Wallraff A 2014 Nat. Commun. 5 5186
[34] Hoi I C Kockum A F Tornberg L Pourkabirian A Johansson G Delsing P Wilson C M 2015 Nat. Phys. 11 1045
[35] Inomata K Lin Z Koshino K Oliver W D Tsai J S Yamamoto T Nakamura Y 2016 Nat. Commun. 7 12303
[36] Forn-Díaz P García-Ripoll J J Peropadre B Orgiazzi J L Yurtalan M A Belyansky R Wilson C M Lupascu A 2017 Nat. Phys. 13 39
[37] Liu Y Houck A A 2017 Nat. Phys. 13 48